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Abstract 

Late time cosmological solutions for scalar field model with kinetic and Gauss 
Bonnet couplings are considered. The quintom scenario is realized with and 
without Big Rip singularity. We find that under specific choice of the Gauss 
Bonnet coupling, the model considerable simplifies, giving rise to solutions 
where the kinetic term is proportional to the square of the Hubble parame- 
ter. This allows to reconstruct the model for a suitable cosmological evolution. 
We considered a solution that matches the observed behavior of the equation 
of state, while Big Rip singularity may be present or absent, depending on the 
parameters of the solution. Evolutionary scenarios known as Little Rip, have 
also been considered. 
PACS 98.80.-k, 95.36+x, 04.50.kd 

1 Introduction 

Recent observational data provide strong evidence for an accelerated expanding uni- 
verse at the current epoch [1], |2], [3], [5]. These astrophysical observations also 
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suggest that the universe is spatially flat with high accuracy, and is composed of 
about 27% of dark matter (including the usual baryonic matter), and 73% of homo- 
geneously distributed new type of negative pressure matter, called dark energy, which 
causes the current accelerated expansion. Among the variety of dynamical approaches 
to the dark energy (DE) problem (see [6], [7], [8]), the scalar tensor theories which 
interpolate between the scalar field and modified gravity models, allow in principle 
to evade the coincidence problem, and in some cases give rise to quintom behavior 
[9], [in]. In the present work, we consider a string and higher-dimensional gravity 
inspired scalar field model, with two kind of couplings: kinetic coupling to curvatures 
and Gauss Bonnet coupling to the scalar field. The GB term is topologically invariant 
in four dimensions and does not contribute to the equations of motion. Nevertheless 
it affects cosmological dynamics when it is coupled to a dynamically evolving scalar 
field through arbitrary function of the field, giving rise to second order differential 
equations of motion, preserving the theory ghost free. Therefore, the coupled GB 
term seems as a natural generalization of the scalar field with non-minimal kinetic 
coupling to curvature [HI |T2] . 

Some late time cosmological aspects of scalar field model with derivative couplings 
to curvature have been considered in [13], [TH [T21 Hlj, [U]- On the other hand, the 
GB invariant coupled to scalar field have been extensively studied. In [16] the GB 
correction was proposed to study the dynamics of dark energy, where it was found 
that quintessence or phantom phase may occur in the late time universe. Accelerating 
cosmologies with GB correction in four and higher dimensions have been discussed in 
[T7] . [T8] . [T9] . The modified GB theory applied to dark energy have been suggested 
in [20], and different aspects of the modified GB model applied to late time acceler- 
ation, have been considered among others, in [21], [22], [23], [21], [2S]- 
All these studies demonstrate that it is quite plausible that the scalar-tensor couplings 
predicted by the fundamental theory may become important at current, low-curvature 
universe (see [26] and [2Z] for review). 

The most general second-order ghost-free scalar-tensor Lagrangian with couplings to 
curvature, can be originated from toroidal compactification of 4 + dimensional La- 
grangian of pure gravity, as shown in [28]. In this compactification the scalar field 
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plays the role of the overall size of the A^-torus, and the couplings appear as exponen- 
tials of the scalar field. This general Lagrangian also appears in the next to leading 
order corrections in the a' expansion of the string theory [29], ^U\. In the present 
study we consider those terms in this ghost-free scalar-tensor Lagrangian that are cou- 
pled to the curvature, where we use more general couplings in order to reconstruct 
appropriate cosmological scenarios satisfying current astrophysical observations. A 
solution unifying early time decelerated behavior, with late time accelerated and 
phantom phases will be studied. We considered solutions with and without Big Rip 
singularity. The equation of state presents minimum at the future, bellow the phan- 
tom divide, turning back asymptotically to de Sitter phase. New solutions known as 
Little Rip, have also been considered. In section II we introduce the model and give 
the general equations, which are then expanded on the FRW metric. In section III 
and IV we present solutions for the Hubble parameter with Big Rip and Little Rip 
singularities, and reconstruct the model according to these cosmological scenarios. 
Concluding remarks are given in section V. 



2 Field Equations 

Let us start with the action for scalar field with kinetic terms non-minimally coupled 
to curvature and coupled to Gauss Bonnet (GB) curvature 

S = [ d^xy/^l-^R - ld^cl)d^^ + F^{<P)G^,d^^d''<i) - V{<i)) + F2((/))d (2.1) 
J LiovrG z J 

where G^i, = R^i, — ^g^^R, Q is the 4-dimensional GB invariant Q = R^ — AR^^R^" + 
RiiupaR^^^" ■ The coupling -Fi(0) has dimension of {length^, and the coupling -^2(0) 
is dimensionless. Besides the couplings of curvatures with kinetic terms, one may 
expect that the presence of GB coupling term may be relevant for the explanation 
of dark energy phenomena. The GB coupling has the advantage that does not make 
contributions higher than second order (in the metric) to the equations of motion, 
and therefore does not introduce ghost terms into the theory. Hence, the equations 
derived from this action contain only second derivatives of the metric and the scalar 
field, avoiding problems with higher order derivatives 
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Taking the variation of action (2.1) with respect to the metric, we obtain a general 
expression of the form 

Rfiu ~^g^jLvR — (2-2) 

where = SttG, is the usual energy-momentum tensor for matter component, 
the tensor T^^ represents the variation of the terms which depend on the scalar field 
and can be written as 

T,, = T^t + Tl^u + T^J" (2.3) 
where T^^, correspond to the variations of the standard minimally coupled terms, 
comes from the kinetic coupling, and comes from the variation of the coupling 
with GB. Due to the kinetic coupling with curvature and the GB coupling, the quan- 
tities derived from this energy-momentum tensors will be considered as effective ones. 
The variations are given by [311 132] 

= V,ct>VA - ^^7^.Va0V> - g,,V{ct>) (2.4) 
= [r,. - \g,uR^ Fi(0)Va0VV + ^7;..VaV^ (Fi(0)V,0V^0) 



- ^(V^V. + V.V^) (Fi(0)Va0V^0) + i?Fi(0)V^0V, 



(2.5) 



- 2Fi(0) (i?,.AVVV.0 + i?,AVVV;.0) + g^,Rx^F^{<p)y^<pV^<p 
+ VaV^ (Fi(0)VVV.0) + VaV, (Fi(0)V"0V^0) 

- VaV^ (Fi(0)V^0V.0) - g^,VxV^ (Fi(0)V>V» 

and 

=4([V^V,F2(0)]/? - g^A^pVPF^mR - 2[W ^F^^R^p - 2[VpV ^F^^R^p 
+ 2[VpVPF2{(t>)]R,u + 2g^,[VPV''F2{ct>)]Rp„ - 2[W^ F^mR.pu.) 

(2.6) 

In this last expression the properties of the 4-dimensional GB invariant have been 
used (see [16], [S])- Variating with respect to the scalar field gives the equation of 
motion 

^ --d^. y^g {RFi{<p)d^ct> - 2R^''F,{cp)dA + 5^0)] + ^+ 



rIF rIF 

— i {Rd^ci)dy - 2R,,dyd''(t)) - -}g = 
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Considering the spatially-flat Friedmann-Robertson- Walker (FRW) metric, 

ds"" = -dt" + aitf {dr^ + r^d^^) (2.8^ 
And assuming an homogeneous time-depending scalar field , the (00) and (11) 



components of the Eq. (2.2), from (2.3 2.6) take the form (with the Hubble parameter 
H = a/a) 



K K 



+ V{(l)) + 9H^Fi{(P)^^ - 2AH- 



■ dFo 



(2.9) 



and 



2H 2Fi 



+ 



+ 8H 



V{(j)) - (3H^ + 2H^ Fi 

• dFy ■; ^ ^^od'^ F2 '2 1 TT TT^^'i 



9 •■ o ' 



+ 16M- 



■ dF2^ 

d(f) . 
(2.10) 

I the present study we have assumed scalar field dominance (the matter term is absent 
in the action). The equation of motion for the scalar field ( 2.7[ ) takes the form 

dFi -2' 



• dV , , 
+ 3H4> + — + 3H^{ 2Fi 



I2HHF1 



24 (^HH^ + 



+ 
dF2 



+ 18H'FA4>U+ 



(2.11) 







where the first three terms correspond to the minimally coupled field. 

Assuming an asymptotic behavior of the scalar field as = 00 = const., and -^2(0) = 

const., then independently of Fi{(j)) he model presents de Sitter solution, as can be 



seen from Eqs. (2.9) and (2.11). From (2.9) and (2.11) it follows that V = Vq = const 
and H = Hq = ka/Vq/S. 



In the following sections we study cosmological solutions of Eqs. (2.9) and (2.11) 
giving rise to accelerated expansion, including quintom behavior, and presenting Big 
Rip and Little Rip singularities. 



3 Solution with and without Big Rip singularity 



We continue studying some solutions to Eqs. (2.9) and (2.11), in the important case 



when the scalar field potential is absent (i.e. V = 0), which leaves the two couplings 
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Fi{(f)) and -^2(0) as the degrees of freedom that describe the cosmological dynamics. 
Making 0^ = tp, from Eq. (2.9) with \^ = 0, it follows 

1 



r-iW = H — H 

3 18H^ 3 dt 



{3.1] 



replacing in (2.11) after simplifications, we obtain 



dH 



>dH dF2 



+ GH^t/j -tlJ^ + 48H^- - 

dt dt dt dt 



dt 



+ 2AH' 



d^F2 .„„odH 



dt^ 



+ 12H' 



dt 



+ 18H^ = 



here we have set = 1. Let's consider the GB coupling as proposed in [32] 

dF2it) g 



dt 



H{t) 



(3.2) 



(3.3) 



Replacing (3.3) in (3.2) one obtains 



dih r, dH , , c,dH , , . 

+ - ^— + 12 2^ + l)H^— + 18 % + l)H^ = 

dt dt dt 



(3.4) 



An interesting and viable late time cosmological solution to this equation is the fol- 
lowing 

(3.5) 



t 7t + ?7 

where p, /3, 77 are dimensionless constants, and a, 7 are constants with inverse time 
dimension. At early times the first term dominates, giving rise to power-law behav- 
ior H p/t. At late times the second term offers an interesting alternative to the 
cosmological constant, as it behaves as a/7, which is approximately the current cos- 
mological behavior. Depending on the sign of 7 and r] may present Big Rip singularity. 
Integrating with respect to time, one finds the scale parameter as 



a{t) = tPe"*/^(7t + r]) 



13 ar] 



(3.6) 



Note that a = at t = 0, indicating that the solution can cover the infiationary 



epoch. Replacing in (3.4) one obtains 



t '-ft + rj 



(3.7) 
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where we have simphfied the solution by setting the integration constant equal to 
zero. Note that in this specific case ?/'(t) = 6H{ty. Integrating the square root of ip 
we find the scalar field 



0(t) 



V 7 
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Note that if the parameters in (3.5) are related by /3'y — ar] = 0, then the scalar field 



simplifies and the Hubble function takes the form H = p/t + P/rj, which has been 
already considered The GB coupling is obtained from ( 3.3[ ) 

'2at + -fp + /3'' 



7t {^A — ajSr] — Spa^rj) 

\ o /A =^ arctan 

a a^'y/Apar] — A 



8a^ 



PI 



y/Apar] — A 
2 /37) In {af + {(3+ p-f)t + pr]) 



(3.9) 



where A = (3 + 2p(3'j + p 7 . And the kinetic coupling from (3.1 ) is 

-2 



Flit) = 



p ^ at + (3 
t 7t + ?7 



(3.10) 



The effective EoS is obtained from the solution (3.5) 



2[(p7^ + /37 — a'r])t'^ + 2p^T]t + pr]"^ 



3{at^ + (p-f + P)t + pr]y 



This EoS has two important limits: 



lim w, 



2 

3p 



(3.11) 



(3.12) 



the particular value p = 2/3 gives the known matter dominance regime. At late times, 



according to Eq. (3.5) the universe becomes dominated by the cosmological constant 
and ends (t — )■ 00) in a de Sitter phase w^fj = —1. We can use the parameters 



appearing in (3.5) to meet the behavior demanded by the current astrophysical ob- 



servations. Normalizing the cosmological time so that the current time corresponds 
to t = to = 1, and taking the current value of the EoS We//(1) = wo = —1 (which 
indeed, according to observations could be a very close to —1) we find the following 
restriction on the parameters 



Pil + + l3'y - a7] = Q 



(3.13) 
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we can also normalize the scale parameter a{t) so that a(l) = 1. This gives the 
restriction 



/3 



a 



(3.14) 



^7 ln(?7 + 7)_ 

these two restrictions allow to express the EoS in terms of the three parameters p, rj 
and 7 as follows 

274r/((27 + r/)t + 1) 



eff 



-1 



(3.15) 



3p [(7 + ?7)2(7t + r])t log(7 + 77) + 777(7 - (27 + r])t)f 

The transition from decelerated to accelerated phase takes place when Weffittr) = 
— 1/3, which gives the following expression to find ttr for given p, 7 and 1] 

p [r/72 - 7/7(27 + 7])ttr + (7 + Hi + V){lttr + V)ttrY = l^V ((27 + V)ttr + Tj) (l-^ir) 

(3.16) 



where we replaced a and /3 from (3.13) and (3.14). There is one more phase presented 



in this solution, which according to the restriction (3.13) is taking place currently: 
the transition to phantom phase delimited hy wq = —1. Then we can expect that at 
t > 1 the EoS takes values We// < — 1. An important feature of the EoS is that it 
does not decrease monotonically with time, but has a minimum at tmm satisfying 

ln(7 + r,) (7(27 + r/)t^,„ - 37'tL„ - 37^^™^ - r/^) + 777 = (3.17) 

As we will see in a numerical example, this minimum occurs bellow the phantom 
barrier, but then the EoS starts increasing towards the value = — 1 at t — )■ 00. 



In the solution (3.5) we can distinguish two important cases: if we consider that 7 



and 7] are positive, then the cosmological evolution avoid Big Rip singularities, and 
on the contrary, if we consider 7 < 0, then the late time evolution faces Big Rip 



singularity. Taking for instance p = 2/3, and using (3.13) and (3.14), we define the 
behavior of the EoS so that at t = 1 crosses the phantom barrier, where 7 and rj are 
used to properly determine the time of transition deceleration-acceleration, according 



to Eq. (3.15). Taking for instance the values 7 = 0.93 and r] = 0.8, then using Eqs. 
(3.13)-(3.17) we find ttr ~ 0.64 and tmm ~ 1-5. In Fig. 1 we show the behavior of the 



EoS in this case 
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0.0 



Fig. 1 The effective EoS for p = 2/3, '-f = 0.93, r/ = 0.8. Note the presence of the 
three phases: matter dominance (decelerated phase) above the line w^ff = —1/3, the 
accelerated-quintessence phase between the lines We/f = — 1/3 and Weff = — 1, and 



phantom phase bellow the line w, 



eff 



-1. 



An important aspect of this solution is that despite the fact that the universe enters 
in the phantom phase, nevertheless the effective energy density remains finite, as can 



be seen by the behavior of the Hubble function (3.5). Similar behavior for the EoS 



was obtained in a quintom model with spinor field [36]. 
Big Rip Singularity 



According to solution (3.5) the future Big Rip singularity occurs provided that 7 < 0. 
Thus, we define ts = ^is the time when the Big rip singularity is approached. 

From (3.5) it follows that as t ^ tg, H ^ 00 and H — )■ 00, meaning that pe// — ^ 00, 
Peff —7- 00, and from (3.6) a — ?■ 00. In Fig. 2 we show the EoS for different cases with 



Big Rip singularity, showing different behaviors of the EoS depending on 7 and rj 




(7;,7)=(l..-0.66) 
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Fig. 2 The effective EoS for p = 2/3 and different values 0/ (r^, 7) .-(1.2, —0.6), 
(1, —0.65), (1, —0.66). Despite the Big Rip singularity, the EoS has a smooth 
behavior. The curve (1, —0.66) returns back to the decelerated phase, entering again 
in the accelerated- quintessence phase at far future. Note that the Big Rip singularity 
may be delayed by making tg = —rj/'-f bigger. 



Therefore, the Hubble parameter given by (3.5) gives rise to Big Rip (BR) singularity 



(Fig. 2), depending on the values of 7 and rj. Thus if rj = —37, then the BR 
singularity occurs at tg = 3, i.e. approximately within 28 Gys. from now. 
Is worth to note that the only future singularity in a, H and H occurs simultaneously 
at t — )■ = —vh (with 7 < 0), when p — )• 00 and \p\ — )■ 00. This excludes type III 
singularity according to the classification given in ([37]), as it requires finite 7^ at 
t — )■ ts- The type II singularity requires finite a and p, while |j9| 00, which can not 
be achieved in the present solution. On the other hand, by simultaneously solving 
the equations H = and H = (keeping 7 < 0) we can find a finite time tc (in terms 
of the parameters a,/3,7 and 77) when p — )■ and \p\ — )■ 0, and the EoS diverges. 
Nevertheless, at this time {tc) the higher order derivatives of H become finite, and 
therefore there is not singularity of the type IV as described in 



4 Little Rip solutions 

In the new solutions known as Little Rip, neither the scale factor nor the energy 
density become infinity in finite time. As in the BR singularity, such solutions cause 
the effect of the structure disintegration in finite time, which can be either earlier 
or later than in a BR model [38], [39], [IQ], [Hj. Let's consider the behavior of the 
Hubble parameter 

Hit) = Hoe^' (4.1) 

where h is a positive constant. Integrating this equation we find the scale parameter 
as 

a{t) = aoe"" (4.2) 
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which leads to the behavior of the density as p oc (Ina)^, which increases with a. 



From (4.1) it follows that H > 0, which reproduces a super-accelerated phase free of 



future singularity. The properties of dissolution of structure of this kind of solutions 
is given in [38] . 

In the frame of the present model we may reconstruct the couplings and the scalar 



field, according to the solution (3.17). Replacing this solution in (3.4) one obtains 
(for the case of 5^ = — 3/4) 



ilj{t) = mle^^' + Ce 



■+ht 



where C is the integration constant. The scalar field is given by 



0(t)= / {%Hle^^' + Ce-"-^+^'Y'^dt 



(4.3) 



(4.4) 



Setting C = 0, the scalar field takes the simple form = y/QHoe^^ /h. The GB 



coupling, as follows from (3.3) is given by 



-ht 



AHnh 



The kinetic coupling from (3.1), and in the case C = 0, becomes 

1 



-2ht 



3Hi 



(4.5) 



(4.6) 



In terms of the scalar field the kinetic and GB couplings are given by the expressions 



3^6 



(4.7) 



^2^ ' 4/l2^ 

Recovering the Newtonian coupling (k^ = SttG), we find that the reconstructed model 
takes the form 



S 



(4.8) 



According to [3HIES], the LR solution (4.1) resembles the ACDM at low redshift and 
is consistent with current supernova observations (which constrain the parameter h). 
In fact we can give a dynamical interpretation to a class of solutions that satisfy the 
requirement of Big Rip or Little Rip (see [3H]-[in]), in the frame of the present model 
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with kinetic and GB couplings. First note that the Eq. (3.4) for the fixed g = —3/4 
takes the form 

+ fei/^ - ^ V - 6^^^^ - 36if^ = (4.9) 
dt \ dt J dt ^ ^ 

then, for any given Hubble parameter H{t), the particular solution 



(4.10) 



satisfies this equation automatically. The corresponding scalar field is found as 



(j)(t) = y H{t)dt 



(4.11) 



The GB coupling is found by integrating Eq. (3.3) and the kinetic coupling from 

1 



(3.1 ) becomes 



Flit) 



3H{ty 



(4.12) 



therefore, for this particular choice of the GB coupling given by (3.3) with g = —3/4 



we can give a dynamical interpretation to any given in advance Hubble parameter, in 
the frame of the present scalar field model with kinetic and GB couplings. In fact a 
Little Rip solution is characterized by a non singular energy density p which increases 
with the scale factor a. Thus, for example taking the scale factor of the form [38j 



a{t) 



(4.13) 



where a{t) is a non singular function of time, the Hubble parameter takes the form 

• 2 

H{t) = a{t). Then, the energy density becomes p = 3a{t) , and the condition 
for increasing density is H = a > 0, which automatically gives super-accelerating 



behavior. From Eq. (4.9 4.11) follows the expression for scalar field 

(f){t) = \/6a{t) 



The GB coupling is found as 



dt 

ait) 



(4.14) 



(4.15) 
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and the kinetic coupling from (3.1) becomes 



3a(ty 



(4.16) 



Of course, from this class of LR solutions we should pick up those that satisfy the re- 
strictions imposed by astrophysical observational data. Particularly at current epoch, 
such models should give a behavior close to the ACDM. 



The reconstruction considered here is based on the particular solution of the Eq. (4.9 ) 



given by the Eq. (4.10). Thus, for a suitable given a{t) we find the scalar field and 



the couplings Fi and F2 through Eqs. (4.14 4.16). Note that this reconstruction is 



applied to late time cosmological evolution. Let us consider the following late time 
dependence of the scale factor 



a{t) 



.At" 



where A is a positive parameter and n > 1. Then, a{t) = At", H 



a 



VeAt". From ( |4.15|4.16| ) it follows 

3 

-t 



2-n 



1 



-t 



2-2n 



' AXn{n-2) ' ' 3A%2 
expressed in terms of the scalar field, these couplings take the form 



AXn{n-2) \^/QX 



(2-n)/n 



1 



3A%2 



2(l-n)/n 



(4.17) 

n\f~^ and 

(4.18) 
(4.19) 



this completes the reconstruction. From (4.17) and the expression for the Hubble 



parameter it follows that the effective density behaves as p oc (Ina)^'^^"^/"'^ which 
for n > 1 avoids the BR singularity and satisfies the condition for LR [38], |39] (i.e. 
ci = if > 0, producing super-acceleration, but the universe needs an infinite time to 
reach p — )■ 00). Note that n = 1 gives the de Sitter universe {p = const.), and when 



n » 1 the couplings behave as Fi 
(ra. 



and F2 



like the couplings given in 



5 Discussion 

In the frame of the scalar field model with kinetic and Gauss Bonnet couplings, we 
considered solutions to the dark energy problem that describe a quintom behavior, 
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with and without BG singularities. According to this solution, it is possible to enter 
into the phantom phase without facing future BR singularities as the case described 
in Fig. 1; but there is also the possibility of getting into a BR singularity, depending 



on the values taken by the parameters, as shown in Fig. 2. The new solution (3.5) 
for the Hubble parameter evolves through the three phases: decelerated expansion, 
accelerated expansion and the final phantom phase with super accelerated expan- 
sion. The time of transitions can be set by adjusting the parameters of the solution. 
Particularly in the examples we considered here, at the current epoch the Universe 
is undergoing the transition to the phantom phase. The results presented in Figs. 
1 and 2 indicate that knowledge of the asymptotic behavior of w{t) (at t — )■ oo) is 
insufficient to distinguish models with a rip from models which are asymptotically de 
Sitter. We have also studied another type of solutions called Little Rip, that cause 
the dissociation of structures in the Universe, due to the increasing behavior of the 
dark energy density. This is the same effect caused by Big Rip solutions, but without 
finite time singularities. 

A remarkable aspect of the present model is that under the choice of the GB coupling 
dF2/dt = g/H with g = —3/4, the equation of motion is automatically satisfied for 
the particular form of the kinetic term as ip = 6H^. For this particular solution, 
the Friedmann equation gives the kinetic coupling of the form Fi = — l/(3if^). This 
allows to reconstruct a class of cosmological scenarios explaining dark energy in the 
frame of the present model, as the examples considered above. 
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